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Abstract. 

We consider generating functionals for computing correlators in quantum field theories with 
random potentials. Examples of such theories include cosmological systems in context of the string 
theory landscape (e.g. cosmic inflation) or condensed matter systems with quenched disorder (e.g. 
spin glass). We use the so-called replica trick to define two different generating functionals for calcu¬ 
lating correlators of the quantum fields averaged over a given distribution of random potentials. The 
first generating functional is appropriate for calculating averaged (in-out) amplitudes and involves a 
single replica of fields, but the replica limit is taken to an (unphysical) negative one number of fields 
outside of the path integral. When the number of replicas is doubled the generating functional can 
also be used for calculating averaged probabilities (squared amplitudes) using the in-in construction. 
The second generating functional involves an infinite number of replicas, but can be used for calcu¬ 
lating both in-out and in-in correlators and the replica limits are taken to only a zero number of fields. 
We discuss the formalism in details for a single real scalar field, but the generalization to more fields 
or to different types of fields is straightforward. We work out three examples: one where the mass 
of scalar field is treated as a random variable and two where the functional form of interactions is 
random, one described by a Gaussian random field and the other by a Euclidean action in the field 
configuration space. 
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1 Introduction 

The standard assumption in quantum field theories is that the quantum action is some fixed funcfional 
of quanfum fields. Similarly in sfafisfical mechanics if is usually assumed fhaf fhe Hamiltonian is 
some fixed funcfional of random fields and ifs conjugafe momenfa. In bofh cases fhe randomness 
in experimenfal oufcomes may only come from fhe fad fhaf fhe (quanfum and/or random) fields are 
described by distributions encoded entirely in the (quantum and/or statistical) partition functionals, 
but there is no randomness in the theory itself (e.g. Lagrangian or Hamiltonian densities are fixed). 
This is, perhaps, a good approach if one is allowed fo perform (and compare againsf each ofher) 
a large number of experimenfs so fhaf fhe inifial randomness in a fheory (due fo fhe fad fhaf fhe 
fheory is nof known a priory) becomes negligible.^ There are, however, af leasf fwo cases when 
fhe uncerfainly in fhe fheory cannof be reduced beyond a cerfain limif and fhus fhe assumption of 
a fixed fheory should nof be used. The firsl case is when fhe fheory confains too many paramefers 
(describing for example a very complicafed potential) all of which cannof be fixed simulfaneously 
due fo experimenfal limifafions. The second case is when fhere are too few experimenfs fhaf one can 
perform on fhe system (due fo, for example, cosmic variance). The firsl case is usually sludied in 
fhe confexl of disordered syslems in condensed matter and fhe second case is presenf in cosmology, 
since fhere is only one universe fo observe. For example, if was known for some lime fhaf a particular 
dislribulion of random polenfials gives rise fo fhe phenomena of Anderson localization [1-4], buf 
more recenlly fhere has also been a quesl fo undersland quanfum inllalionary perlurbafions on random 
potentials [5-7]. The silualion is complicafed further in fhe confexl of siring fheory models of elernal 
infialion, where one has to deal wilh bofh cases simulfaneously: a large number of paramefers (of 
fhe siring fheory landscape [8-10]) and only a single observational evenl of fhe Universe we live 
in.^ To delermine fhe sfafisfical properfies in such syslems one usually generales numerically a single 

^Note, that the initial uncertainties would never get completely eliminated, but they can be made arbitrary small with 
more and more data coming from experiments. 

^Strictly speaking the models of eternal inflation also suffer from an additional uncertainty of choosing a probability 
measure [11] (known as the measure problem) which will be discussed in a follow-up publication [12]. 
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random potential and then simulates directly the classical background dynamics [13-19]. In this work 
we present a fully analytical (and also quantum) treatment of random potentials using the framework 
of path integrals by defining two generating functionals both of which can in principle be used for 
calculating averaged correlators. The analysis is done for a single real scalar quantum field theory 
with a distribution in potentials, but the generalization to many fields and to different types of fields 
is sfraightforward. 

By random potential we mean that the potential had been drawn from a statistical ensemble 
described by some probability distribution P[T] which is known, but the exact realization of the po¬ 
tential V (tp) is unknown. The simplest case perhaps is if there are no (non-quadratic) interactions, i.e. 
V (tp) = and all of the randomness is imprinted in the value of the mass m, chosen according 

to some distribution P(m). A distribution in the coupling constant of the non-quadratic interactions 
may be chosen as well. A more interesting case though might be to consider a functional distribution 
for the potential which can be described by some (Euclidean) action in field’s configuration space. 



with 



(1.2) 


where hy is an arbitrary constant in the distribution P[y] that need not be related to Planck’s constant 
h=l. 

Note that the potential ‘field’ V (tp) lives in the configuration field space and not in physical 
space-time and the expression of Eqs. (1.1) and (1.2) is nothing but a Euclidean path integral of 
‘field’ V with ‘time’ variable tp. A discretized version of Eq. (1.1) is understood as 


where s is the discretization scale of ‘time’ variable tp, V (tp,) are the values of ‘field’ at discrete points 
tp,, and the measure of path integration in (1.2) is understood as 

DV = Y[dVi. (1.4) 


This in turn will make the integral in the path integral over field (i.e. /D(p/[(p] = /n;JtP//(tP 0 )tPi 5 •■•)) 
to be rather a sum 

[ W[(P] (1-5) 

'0 h h 

where /[tp] is an arbitrary functional of function of space-time tp(x). This shows how the discretiza¬ 
tion can be carried out in principle, but we will attempt to work in a continuum limit for the remainder 
of the paper. 

As discussed in Sec. 3.1, having ‘Eagrangians’ for potential described by Eqs. (1.1) and (1.2) 
leads to a technical issue and analytical calculations seems formidable, but one can still tackle the 
problem numerically. More generally, the statistical ensemble of random potentials need not be 
described by a local action as in (1.1), for instance it could be described by some Gaussian random 
field with a prescribed two point correlator Ay for the potential, i.e. 


P[V] = ;^exp (^-JJ r/(pir/(p2y((pi)A^,'((pi,(p2)y((p2)j , 


(1.6) 
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where the integration can be understood by discretizing the field space as was described in the previ¬ 
ous paragraph if needed. In either case, the methods of generating functionals developed in this paper 
apply. 

The properties of quantum (or random) fields such as amplitudes or correlators are conveniently 
encoded in generating functionals for correlators also known as quantum (or statistical) partition 
functional. For example, the partition functional for a real quantum scalar field is given by [20, 21] 



Dtp exp 




(1.7) 


where we have set = 1. Then the in-out ^ correlators can be obtained by functional differentiation 
with respect to sources at different points, with proper normalization, i.e. 


(r{(p(xi)..(p(xi)}) 


{-if5^Z[V,J] 


Z[F,0]67(x^)..57(xi) 


7=0 


( 1 . 8 ) 


where T stands for time ordering. From now on, we suppress writing out the T operation, but it is 
always assumed to be present inside correlators. Next, the generating functional for in-in correlators ^ 
is obtained by adding one more (conjugated) partition function (i.e. Z[F,7+]Z* and thus having 

two fields (p+ and tp^ with respective sources 7+ and J^), and the in-in correlators can be obtained 
by functional differentiation with respect to sources at different points, i.e. 




5*^z[F,7+] 

Z* [F, J ]5J (y 1)..67- (yt) Z[F,7+]57+ (x^) ..67+ (xj) 

(1.9) 


Note that in (1.7) it is assumed that F((p) is a known potential. The question that we would like to 
ask in this paper is what should be the generating functional if the exact form of V (tp) is unknown or, 
in other words, if the masses (or any other parameter in the Lagrangian) and/or the functional form 
of the potential/interactions are random. Note that the same question can be asked in the context of 
statistical physics systems where this additional randomness may come from the Hamiltonian itself 
[3, 4, 20, 27]. 

There are four types of averages that we will have to distinguish. Two being statistical in origin 
(known as annealed and quenched averages in statistical physics [3, 4]) and two being quantum in 
origin (in-out and in-in averages/correlators in quantum field theories as recapitulated above). The 
statistical averages will be denoted with a bar (Oa for annealed and Oq for quenched) and the quantum 
averages will be denoted with square brackets {{O) for in-out and {0)*{0) for in-in) as usual, where O 
is some operator (observable). The distinction between annealed and quenched averaging or between 
in-out or in-in averages should be clear from context. 

Annealed average must be taken whenever the target field tp can backreact on the random field 
V and both evolve together (for example as in regular interacting field theories with V replaced by 
another field ’living’ in the same physical space-time).^ Then the averaged in-out correlators, for 


^For most of this paper, we refer to in and out vacuum states (i.e. \vac,in) and \vac,out)) as just in and out states 
respectively without mentioning the term vacuum, with the hope that it is clear from the context. 

^We refer to in-in as multiplied in-out and out-in correlators/amplitudes, where out-in are nothing but conjugated in-out 
correlators, i.e. (in, vac|r^ {(p...(p}\vac,out) = (ok?, vac|r{(p...(p}|vac, in)*. Although these are not exactly the standard in- 
in objects as usually defined in QFT’s [22-24] and cosmology [25, 26] where a sum over all possible out states is implied, 
we will still use the term in-in with the understanding that there is only one out state (the vacuum out state). 

^There are many other places where annealed averages are studied, e.g. spin glasses in condensed matter physics[3, 4, 
27, 28]; quantum brownian motion[20, 29, 30], also see [31]; and non-equilibrium QFT[20, 32-34]. 
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instance, are given by 


((p(v^)...(p(vi))^ 


/DVP[V]/D(pc'Xvi)...(pfa) 


( 1 . 10 ) 


where k denotes what point correlator it is. Note that since both the random potential V and the field 
tp can be dynamic, we simply inserted one more path integral for the random field V. Similarly, the 
corresponding averaged in-in correlators (1.9) are 


((p(y^)...(p(yi))*((p(;ci:)...(p(;ci))„ 


jDVP[v]fD(?e X(yj^)-tp(yi)/f^XX-^i)-tpfa) 
/ DVP[V] f D(pe-‘^ f D(pP^ 


On the other hand in the case of quenched systems, the random field is quenched (or frozen) and does 
not evolve in time. This means that the target field whose evolution depends on the random field does 
not backreact on the latter, which is the case we will be addressing in this paper. Since our random 
field manifests itself as a potential V (tp) in the Lagrangian for the target field tp, we want to study 
field theory on random, but with a given distribution P[V], potentials in order to obtain the averaged 
correlators. Hence, the quenched averaged in-out correlators are given by 


«»)>, = JdvpIv] 


( 1 . 12 ) 


and correspondingly the quenched averaged in-in correlators are 


(1.13) 

It is understood that the action S for the field tp depends on the potential and therefore we need not 
bother writing it as a functional of V over and over again. Since in this paper we are interested in 
quenched averages, we do not use annealed averages from here on, moreover O without a subscript 
always indicates a quenched average. 

This paper is organized as follows. In Sec. 2 we summarize the replica trick and define a 
generating functional for quenched averages with a single (or double) replica of fields for calculafing 
in-out (or in-in) correlators. In Sec. 3 we use the generating functional for calculating averaged 
correlators for a distribution of random self-interactions described by either a partition function of a 
Euclidian harmonic oscillator (see Sec. 3.1) or by two-point correlators of a Gaussian random field 
(see Sec. 3.2). In Sec. 4 we construcl a generating functional for quenched averages with an infinite 
number of replica of fields which we use in Sec. 5 fo calculate averaged correlafors for a particular 
distribution of random masses. The main results of the paper are summarized and discussed in Sec. 
6 . 


2 Single (Double) Replica(s) of Fields 

The main challenge in calculating quenched averages is the dependence on random potentials in the 
denominators in expressions (1.12) and (1.13). The problem is usually handled by either the so-called 
replica trick or by introducing Grassmanian fields, see e.g. the study of disordered systems in con¬ 
densed matter [2-4, 20, 27]. However one major difference from those condensed matter studies is 
that there the potential ‘lives’ in the same physical space where the field ‘lives’, whereas in our case, 
which is relevant for cosmological systems and in general for relativistic systems, the potential ‘lives’ 
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in the field space of tp. Also, we want to obtain a generating functional to generate all averaged corre¬ 
lators, all of which are potentially observable quantities. In what follows we present two generating 
functionals for generating averaged correlators, with one (see Sec. 4) being superior and more trust¬ 
worthy compared to another in the sense that only interpolation of the functional between positive 
integers is needed and not extrapolation to negative integers (discussed later), but at the expense of 
increased complexity. The formalism is illustrated for a real scalar field with random potentials (see 
Sec. 3), and random masses (see Sec. 5), and we note that the generalization to different types and 
collections of fields is straightforward. Possible application of the proposed recipe to inflationary 
cosmology will be discussed in detail in a separate publication [12]. 

To illustrate the replica trick we can make use of a trivial identity 

zjibi 

where Z[y,7] is the partition functional with a fixed potential. Then generating functionals for quan¬ 
tum systems with random potentials with single and double replica can be defined respectively as 


[/, f\ = lim [ DVP[V]Z [y, J] Z" [y, f\ 


( 2 . 2 ) 


and 

z“[7+,/+,7^,/^] = lim /'Dyp[y]z[y,7+]z”^[y,/+]z*[y,7-]z*"^[y,/-]. (2.3) 

and the quenched averaged in-out and in-in correlators can be obtained by differentiating (2.2) and 
(2.3) with respect to sources, i.e. 


((p(xi)...(p(xjt)) 


(-0^6^ 


dJ{xk)..dJ{x\) 


j=J=o 


(2.4) 


and 


((p(yi)...(p(y^))*((p(xi)...(p(x;t)) = 




57 {yk)-^J (yi)67+(x;t)..67+(xi) 


7±=/±=0 


(2.5) 


^ However, for the last two equations (2.4) and (2.5) to be correct we must make two crucial as¬ 
sumptions. First of all we must assume that the integrands of the path integrals in (2.2) and (2.3) as 
functions of n and ,rr respectively can be not only interpolated between integer values, but also 
extrapolated to negative values. Secondly we must also assume that one can interchange the order 
of taking limits and path integrations in (2.2) and (2.3). In Sec. 4 we will construct an improved 
generating functional W with multiple replicas where the extrapolation assumption can be relaxed 
by avoiding taking limits to negative numbers of fields. 

Moving on, if there is a functional distribution for potentials, e.g. a Lagrangian for self¬ 
interactions (1.1), a closed form integration of the field’s partition function is not possible in general; 
hence, we resort to a perturbative calculation in this case. For this treatment to work, we need to 
assume weak random interactions (i.e. interactions with a small coupling constant g) so that a pertur¬ 
bative treatment makes sense. Further there is no limiting procedure required in (2.2) and (2.2) as n’s 
are equal to negative one with binomial expansion at work. For example, for a real scalar field with 
action 

= j -^m^(^'^ + gV{(^)+J(^ , (2.6) 

®Note that the tilded sources are introduced precisely to handle dependence of V in the denominator, and they are not 
used for functional differentiation to obtain correlators. 
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the partition functional can be expanded perturbatively in g as 


z[y,7] = 



around the free partition function 


Zf[J] = exp 



v)J{v)d^ud^v ) , 


(2.7) 


( 2 . 8 ) 


Here, Af{x — y) is the free field propagator (negative of the inverse of the free field operafor), also 
called Feynman Green’s funclion, and if is equal fo 


Af (x) = lim 


e^O (27I)"^ 


d^k 




— I 


— rrP- + ie (27t)^ 


0(x°)- 






2(0i 


+ 0(-x°)- 




-k.x) 


2a)k 


d^k 


(2.9) 


where fhe energies cOit = + Vk^ + nP- are assumed fo be posifive and and we are using diag(+, —,—,—) 
as fhe Minkowski signafure. The conjugafed propagator (in fhe conjugafed partition funcfion) would 
be, aparf from a negafive sign, reversed in lime: 


A^(x) = 


(271)3 


((cojjrO-f.x) -i((akx°-k.x 

0(x°)^---- + 0(-x°)-^ 


2(0^ 


2(i)k 


d^k 


( 2 . 10 ) 


i.e. fhe posifive and negative frequencies gel swapped wifh respecf fo posifive and negafive times. In 
the following section, we describe how the perturbation theory works for random interactions given 
by some distribution P[y] for a real scalar field wifh fhe action (2.6). 


3 Random Potentials 

As described before, for a perlurbalive Irealmenl in g (i.e. field (|) wifh random self-inferacfions), no 
limiling procedures are required and all n’s in (2.2) and (2.3) can be sel fo negative one. Then fhe 
in-ouf and in-in parfilion funclion can be expanded in g as 

+ (/g)2z“(2)[7,7]... (3.1) 

and 






(3.2) 

where fhe zerolh order ferms are 




II 

Ary Ai 

(3.3) 


Zf[J+]Z}[J-] 

Zf[j+]z}[j-r 

(3.4) 
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and the first order terms are 


= DVP[V] / 


( /D(pV[(p(z)]e'^/M /D(pV[(p(z)]e'^/[^’^']\ Zf[J] 




Zf[J] 


Zf[J\ j Zf[J\ 


(3.5) 


= jnVP^] Jd^z 


f D(p+V[(p+(z)]e'‘^t'['P^.-/+] J D(p+ 


+ 


Zf[J+] Zf[J+] 

J Dtp^ /Dtp^ 2/[7+]Z^[7^] 


Z}[J^ 


Z}[J- 


J Zf[J+]Z}[J-] 


where Sf denotes the free action (i.e. (2.6) with g = 0), and so on. By swapping the space-time 
integral(s) and field pafh infegral(s) wifh pafh infegral(s) in fhe pofenfialj we gef correlators in po¬ 
tential as inferacfions (in such a way fhaf fhe /:-lh order ferm in g has /:-poinl potential correlafors as 
inferacfions). For example, we have for fhe firsf order terms 


^io{l) = fd\( Zy[7] 


and 


/+,/-,/-] = / d\ 


Zf[J] Zf[J\ J Zf[J] 

J D(p+F[(p+(z)]e'‘^t'[9+../+] J D(p+F[(p+(z)]e'‘^t[9+/^] 


(3.7) 


-h 


V Zy-[7+] Zy[7+] 

/D(p-F[(p-(z)]e-'^/[^^7-] /D(p-F[(p-(z)]e-'^/[‘P"7-]\ Z/[7+]Z}[7- 


z;[/-] z;[7-] 

where any potential correlafor is defined as usual 

F((pi)...F((p^) = J DFP[F]F((Pi)...F((pO. 


Z/[7+]Z*[7- 


(3.8) 


(3.9) 


Now fhaf all of the interactions are encoded in correlators of potentials we can carry on the per¬ 
turbative calculations for a given set of such correlators. In what follows, we describe two simple 
examples: in the first one, the desired set of correlators arises when the potential is modeled as a 
random field described by some Lagrangian in fhe field tp configuration space, and in fhe second one, 
we consider a Gaussian random potential described by some fwo poinf correlafor of fhe potential. In 
bofh of our cases fhe one-poinf correlafor vanishes (which mighf nof be frue in general) and fhus we 
need to expand to second order in g in order to see the effect of random interactions. For the in-out 
generating functional the second order term is given by*^ 


’ JJ I 2Zf[J] 2Zf[J\ 

//D(pD(pV[(p(z)]F[(p(z')]e'^^''P’^'+'^^'^’-'’l //D(pD(p'F[(p(z)]F[(p'(z')]e'^^t^’^'l+'^^[^'’^'’'\ Zf[J] 

Zf[J]Zf[J\ + Zf[J\Zf[P] )zf[J\ 

^We assume that for any given potential V (tp) drawn from the distribution P\y\, the space-time integral(s) and field path 
integral(s) converge uniformly and hence the integral(s) can be swapped. 

®Prime is introduced for the squared terms in the perturbation series in the denominator: Since each quantity in the 
perturbation series is an integrated quantity, squaring (or any power in it) must be handled by differentiating between each 
other. 
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and a similar (but much messier) second order term can be derived for the in-in generating functional. 
We first start with the Lagrangian approach where we identify a technical problem which does not 
allow the analysis to be carried out analytically. The problem is due to a discontinues derivative of the 
correlators at zero and is likely to be present for the simplest Lagrangians, but can be easily handled 
numerically. Then we perform a sample analytic calculation for Gaussian random potentials with 
sufficiently smooth two-point correlators. 


3.1 Euclidean Harmonic Oscillator 


The partition functional for potentials (not to be confused with the partition functional for fields (1.7)) 
described by a Lagrangian (1.1) can be obtained by inserting a source Jy coupled to the potential V: 


Zv[V,Jv] = JoVexp Jd(i;> 




(3.11) 


where is the Lagrangian for the potential. In the case of a Euclidean harmonic oscillator we have 


Zy[L,7y] = DVexpl- 


r/tp 


1 fdvy 1 , , 

2 U) 


(3.12) 


Note that a particular limit of this distribution is known to give rise to the phenomenon of Anderson 
localization when hy , my —)• 00 , such that my / hy remains finite. Depending on the values of hy and 
my, we have a Gaussian distribution in both the potential and its derivative (centered at zero) at every 
point in the field space. To carry out the matrix Gaussian integral we make the assumption that 
the potential dies off faster than l/tp on either infinities, so that we can neglect boundary terms. The 
corresponding Green’s function (with a single argument) of the differential operator ^ 
is given by 


my 






dk- 


„ik(f 


This integral can be evaluated without any deformations of contours as the poles {—imy,imy) are 
on the imaginary axes. For tp > 0, we take the upper half plane in the complex k space to close the 
integration contour with the real axis and for tp < 0, we take the lower half plane to get 


Ay((p) = ^ r0((p)c-'”^‘P + 0(-(p)c'"'''Pl = —e 

my ^ ^ my 


(3.14) 


It is straightforward to generate all higher points correlators of the random potential (3.9) from the 
potential’s Green’s function (3.13): all odd point correlators are zero, the two point correlators is 
simply 

E(cpi)L((p 2 ) =Ay((pi-(P 2 ), (3.15) 

and all higher even-point correlators are obtained from (3.15) using Wick’s theorem [20, 21]. 

With that, we can now move back to the generating functionals (3.1) and (3.2). Retaining 
only the leading non-zero correction (which is second order in g for the considered example) and 
replacing the fields as functional derivatives with respect to the respective sources, we get for the 
in-out generating functional (3.1): 




z/ 

{ig? 

2/Z/ 


1 + 


jig? 

IZf 


d^zd^z y 


d^zd'^Z Ay 


f67(z) 

6 


i5J{z'] 


Zf- 


(ig? 

2Zf 


d‘^zd‘^z! Ay 


ihJ{z) ihJ{z') 


2/^/ + 


{igf 

Z/Z' 


d'^zd'^z! Ay 


iBJ{z) 

6 


idJ{z) iBJ'{z') 




ZfZ'f 


(3.16) 






















and a similar expression for the in-in generating functional (3.2). Here we have suppressed the 
source dependence of the free partition functions, hence the tildes Zf = Zf[J\ and primes z/ = Zf[J']. 
The correlator has a derivative discontinuity at zero and therefore replacing the fields as functional 
derivatives w.r.t. sources (as usual in path integral QFT) is problematic at zero. Note that we would 
be facing the same problem for the in-in generating functional (3.2). This can be dealt with by either 
modifying the correlator (smooth it at zero) which would modify the Lagrangian (see next section), 
or by starting with well-behaved correlator(s) which may or may not be described by any Lagrangian. 
Note that there can still be Lagrangians with sufficiently smooth correlators, but it is not the case for 
the example of the Euclidean harmonic oscillator considered here. 


3.2 Gaussian Random Potential 


Consider a Gaussian random potential defined by fhe probabilify disfribufion 


= )^exp 



)AyHcpi,(p2)E((p2) 


(3.17) 


whose odd-poinf correlafors vanish and all even-poinf correlafors can be obfained from fhe fwo-poinf 
correlator Ay using Wick’s fheorem [20, 21]. As an example, we smoofh ouf fhe correlator (3.13) 
using hyperbolic fangenf (i.e. replacing \x\ wifh xlanh(Yx)): 


Ay((p) 


hy 

my 


exp(—|my|(planh(y(p)) = Ay((p) 


(3.18) 


wifh Y being a smoofhing consfanf andwe use fhis expression in fhe Gaussian random disfribufion.^ 
The inferacfions are still described by fhe in-ouf generafing functional (3.16) wifh Ay replaced by 
Ay and we can calculafe Feynman diagrams by expanding fhe smoofhed pofenfial Green’s function 
(3.18) around zero: 

Ay((pi - CP2) = ^ ^1 - my Y(CP1 - ^2)^ + + 1 !}^^ ((p^ _ , (3J9) 

Using (3.19), we approximate (3.16) (wifh smoofhed ouf correlator (3.18)) furfher by fhe leading 
order terms in y which gives rise to quadratic interactions. After all the usual diagram arithmetics 
we arrive at 


2;'° = 


1 + ^-^ II d\d‘^z{ j. -My, 


+ 0(/ 


Zy[0] 


where 


j, _.yy ) '^^^^■7(M)Af-(M — rv). 

with solid dots representing internal vertices. The averaged two point correlator is then 

522'^^ 

7=0 


((P(X2)(P(X1)) = T 


i5J{x2) i5J{xi) 


-i(^xi ^ X2 ) -g'^^vy\ 


z z 

Xi X2 


(3.20) 


(3.21) 


(3.22) 


®Note that this correlator might be describable using the Lagrangian formalism of the previous section (with perhaps 
added correction to the action (3.12) that depends on some negative power of y so that in the limit y —>■ oo, we recover the 
Euclidean harmonic oscillator), but the exact form of the Lagrangian is not important for our purposes. 

*®Since the sourced fields denoted with tildes are never to be differentiated when obtaining averaged field correlators, 
we can set these sources to zero after calculating respective Feynman diagrams. 
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The first diagram is just the usual propagator 


x^y ) ^Apiy-x) 
and the second diagram a correction representing the following integral 


ww' ^ JJ Apiw — x)Af{y— W) 


(3.23) 


(3.24) 


where arrows specify the direction of time flow (i.e. having usual propagators (2.9))*^. Note that 

(3.24) is nothing but a multiplication of two diagrams: creation of a particle at v and its annihilation 
somewhere in between, followed by the creation of a particle somewhere in between and its annihi¬ 
lation at y. The integral can be computed easily since the two diagrams (integrals) are disconnected, 
so that 



(3.25) 


A similar calculation for the in-in generating functional (3.2) gives 


z}' = [i+iM jj 

+ 2( j+, - .21 


7+ 


•Zi 


•Z2 


J+- 

+ 0(/) 


-22 j + 
(Z/[0])2 


•Zl 


‘Z2 


(3.26) 


where the bracketed diagrams with 7+ sources are as described before (eq. (3.21)), and 7 sources 
are the same but with conjugated propagators: 


w ) = 


d^uJ 


{u)A*p{u — w). 


(3.27) 


Thus, for example, the averaged two point in-in correlator reads: 


(cp(y2)(p(yi))*((p(x2)(p(xi)) 




57 (y2)57 (yi)57+(v2)57+(vi) 


j+,j-=o 


j -ig^hvj 

y 2 I 


( 


\ 


Z\Z2 

X\ > • •) X'2, 


\yi 


yi 


Xi ->- X2 

yi y2 

ZlZ2 


(3.28) 


where the time reversed flow (arrows from right to left) specifies conjugated propagators i.e. 

(;c^y)=A*(v-y) (3.29) 

and 

JJ d^wd'^w'A*p{w'— y)A*p{x — w) = — (3.30) 

Also, diagrams appearing one on top of another are nothing but a product of them. Notice that the 
averaged one point in-in correlator is also non zero due to the third term in (3.26) even though all 

"We choose a convention where time flows from left to right. 
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odd-point correlators for potentials vanish, since in totality there are two field operators (causal, anti- 
causal) which is even, i.e. 


{9(j)>'<ipW) 




67 (y) 67+ (x) 


= g 


7+,/-=0 



m‘^ 


(3.31) 


This suggests that both the normal and time reversed fields inferacf wifh each ofher fhrough fhe ran¬ 
dom pofenfial. Since we were doing a perfurbafive expansion no questionable assumpfions (discussed 
above) had fo be made abouf fhe generating functionals (2.4) and (2.5). However for a random dis- 
fribufion of masses P{m) (which we are going fo discuss in Sec. 5) a closed form infegrafion over 
random pofenfials is possible, and one mighf quesfion fhe necessity of extrapolating generating func¬ 
tionals fo fhe unphysical negafive one number of fields. In fhe nexf secfion, we will develop a second 
generating functional which is a lof more symmefric and, moreover, only needs fo be interpolated 
befween (a lof more physical) zero and one number of fields. 


4 Multiple Replicas of Fields 


We sfarf by developing fhe generafing functional for averaged in-ouf correlafors and show how fo 
calculafe fhe corresponding averaged correlators from if. Extension fo fhe in-in generafing functional 
is sfraighfforward. Firsl nofe fhaf fhe issue of having V dependence of partition funclions in denom¬ 
inators of (1.12) and (1.13) can be deal! wifh by using logarifhm which, upon differentiation wifh 
respecf fo fhe sources produces fhe required denominator automatically. Hence we consider fhe aver¬ 
age of logarifhm of fhe partition funcfional (generafing functional for connected Feynman diagrams: 
IT = InZ) 

W[J] = f DVP[V]W[V,J] = lim -^ ( / DTP[T]Z”[y,7] ) . (4.1) 

7 «^o+ dn J 

The exponenfiafed objecf Z" is well defined af all non-negafive integers n, buf fhe gain here is fhaf we 
can sfudy fhe limif of fhe {interpolated) funclion af n —)• 0+(i.e. having no fields). In fhis consfrucfion 
we are avoiding fhe need for exfrapolafion fo negafive values of n which is due fo fhe facf fhaf fhe 
objecf is more symmefric as no disfincfion is made befween partition functions in numerators and 
denominators in (1.12) and (1.13). If (4.1) is functionally differentiated wifh respecf fo sources, fhe 
desired averaged one poinf in-ouf correlator is indeed generated 


5W 

/67(x) 


7=0 


DVP[V] 






(4.2) 


buf fhe generafing funcfional cannof be used to generafe all higher poinf correlafors. Whaf if generates 
fhough are fhe averaged cumulants of in-ouf correlafors, for insfance 


62 IT 


ihJ{x 2 ) ihJ{x]) 


7=0 


((p(x2)tp(xi)) - ((p(x2))((p(xi)) 


(4.3) 


and so on. This is due fo fhe facf fhaf ((p(x 2 ))(tp(xi)) / ((p(x 2 )) ((p(xi)), and fhus we are nof gener¬ 
afing cumulanfs in fhe usual sense. 

We could have been expecfing fhis property, since we are averaging only over fhe connecfed 
diagrams, whereas in general, an in-ouf correlator has bofh connecfed and disconnecfed diagrams. 
Hence, we need a machinery fo build all fhe missing averaged disconnecfed diagrams. To accomplish 
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this task, we propose a generating functional with an infinite number of replicated sets of identical 
fields, each set with a different source*^: 


%^^°{nun2,...)[JuJ2-]= = j DV P^]Z'^^[V,Ji]Z^^^,J2]... (4.4) 


The averaged products of different IT’s are obtained by differentiations with respect to the n’s. 




nin2...nk 


= W[V,Ji]...W[V,Jk] = 


dP 


dni dn2...dnic 


ni 


(4.5) 


where k is an arbitrary integer. These are the averages which will be used for building different 
correlators. For example, to get the averaged one point correlator we only need to differentiate (4.4) 
with respect to any one n and the respective source only once, i.e. 


((p(x)) 


i5Ji (x) 




niJi{x) 


(4.6) 


which is equivalent to (4.2). For the averaged two point correlator, we need two terms (the second 
one is needed to cancel out the second term in (4.3)): 


((p(v2)tp(vi)) = ((p(v2)tp(vi))-((p(v2))(tp(vi)) + ((p(v2))(tp(vi)) = 


(4.7) 




rii 


idJ\ (x2) /57i (xi) 


+ 




«1«2 


7i=0 


ihJ2{x2) fS/i (vi) 


= HA/iO _j_ q^io 

1 /l( j:i) 7l(X2) ' ^'^nin2Ji(xi)J2(x2) 


Jl,J2=0 


and so on. 

There is a pattern that one might notice in calculating all these averaged correlators so that 
we can identify simple rules for obtaining them. To understand the origin for these rules it might be 
helpful to recall from QFT that when the generator for connected diagrams W is functionally differen¬ 
tiated with respect to its source, every of the closed brackets ((p(v/)...(p(vi)) after differentiation with 
respect to source at Xk goes to ((p(vjt)cp(v/)...(p(vi)) — ((p(vjt))(tp(v/)...(p(xi)). Also notice that TF'” 
differentiated with respect to a source without being differentiated with respect to the corresponding 
n and evaluated at zero vanishes. Hence, to get an averaged k point in-out correlator, we need to form 
all partitions of the set containing k space-time points {x\,...,Xk}- The total number of partitions is 
known as Bell’s number; each of the partitions corresponds to a term and all such terms have to be 
added to give the averaged in-out correlator. For example, for k = 3 the total number of partitions 
is five: {{xi,X 2 ,X 3 }), {{xuX 2 },{x 3 }),{{xi,X 3 },{x 2 }),i{x 2 ,X 3 },{xi}) and ({xi},{x 2 },{xa}) and the 
corresponding terms are ‘^m,n 2 ,« 3 ,/i,/ 2 ,P 

tively. It is assumed that the space-time dependence of the sources is ordered such that the first source 
depends on xi, the second source depends on X 2 , etc., e.g. 




H1 H2 7l J 2 J\ 


dnidn 2 5/i (xa) 6/2 (x 2 ) 5/i {x 1 ) 






(4.8) 


By adding all of the terms we get 


((p(xa)(p(x2)(p(xi)) = + 

'^^in 2 JiJiJ 2 + ^^nin2Ji J2J2 n2Ji J2J1 « 2 M 3 .tl J2J2 (4.9) 

*^See [35], section III wherein two replica sets have been used to calculate the averaged product of two one-point 
correlator for statistical systems. 
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So far we have only discussed the one-, two- and three-point correlators, but the same construction 
of averaged correlators from partitions of the set of space-time coordinates {xi,...,Xk} works for an 
arbitrary k. 

Moving on to the calculation of in-in correlators, the generating functional is a straightforward 
extension of (4.4): 


= I DVP[V]Z<[y,7+]...Z*”»'[V,7^]. 


(4.10) 


The two pieces (in-out and out-in) of the averaged in-in correlators in (1.13) are obtained from the 
in-out and out-in integrands of (4.10) respectively as described above. Since each of these two pieces 
are obtained from adding terms corresponding to partitions of the two sets of space-time coordinates 
{x\,...,Xk} and {y\,---,yk} (representing the coordinates of in-out 7+ and out-in J sources respec¬ 
tively), the total number of such terms is equal to the product of two Bell numbers. For example, the 
averaged in-in two point correlator is given by 


((p(y2)cp(yi))*((p(v2)(p(vi)) 




7+7+ 4” 
7l Jy 




i j fi'2^ A/ j j Y j j 


/+ /+ 
^2 Ji 




ft2 At 1 At j ^2 ^ J 


/+ /+ 
Jy Jy 


+ W' 


rt I 




(4.11) 


Since there are two terms for both in-out and out-in individually, there are four terms in total. 

It is worth mentioning that the generating functional for averaged in-out correlators is 
contained within the generating functional for in-in correlators and thus the latter is a more 
general object. To get the former, one should not differentiate with respect to the minus sources and 
follow the regular procedure as described above for the averaged in-out correlators. Furthermore, the 
Euclidean version of this generating functional i.e. one with Euclidean fields, is straightforwardly 
obtained by doing a Wick rotation (t —)• —ix for ’-t’ fields and t —)■ ix for fields). We now move 
on fo real scalar field fheory wifh a random mass and use fhe improved, more general generafing 
funcfional (4.10) fo calculafe averaged in-ouf and in-in correlafors. To illusfrafe fhe idea and be more 
simplisfic, we do fhis for Euclidean scalar field. The resulfs for Eorenfzian/Quanfum field fheory 
can be easily obfained by undoing fhe Wick rofafion in fhe end fo obfain Quanfum correlafors from 
Euclidean correlators (x —>• it and x —)■ —it for ’-t’ and fields respectively). 


5 Random Masses 


To illusfrale fhe effecl of random paramefers in fhe acfion, we consider a massive scalar Euclidean 
field fheory (Wick rofafed Eorenfzian fheory) wifh source 

SE{m)[(^,J] = Jd'^XE Q(d(p)^-h^m^(p2-7(p^ (5.1) 

where 



(5.2) 


wifh fhe following normalized disfribufion for mass 


P{m) = 


1 


a\/^ 


exp 


[m — nio 
2a^ 


(5.3) 
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The subscript E stands for Euclidean which we drop from now on assuming it is understood. Since 
the in-in generating functional (4.10) is a more general object, in the sense that it can generate both 
in-in and in-out correlators, we work with it (Euclidean version offcourse). Eurthermore, to replicate 
all the partition functions in (4.10) we assume that all n’s are positive integers. Accordingly, for a 
scalar field with action (5.1) and distribution (5.3) of masses, the generating functional is given by 



where the second index for the sources (and fields) looks redundant, but is needed, because each of the 
partition functions (within a replicated set) is an integrated quantity and must be handled separately. 
Once we have calculated the partition functions, we drop the second subscript. To keep our notation 
simple, we don’t bother writing = ... = /„ ; = ... = anymore, with the understanding that 

it is implied. 

Since each of the partition functions are divergent quantities if directly calculated on a contin¬ 
uum space, we discretize Euclidean space-time and perform lattice field theory; the properly normal¬ 
ized correlators are obtained in the continuum limit. Coming back to (5.4) with Eucledian space-time 
now discretized (i.e. coordinates x being replaced with labels k of lattice points, e.g. tp^p (x) —)■ (p^p ^), 
we separate the mass dependence from the partition functions and bring the integration over the mass 
inside the scalar fields’ path integrals: 


4E“ = 


where 


a=l \p=l Y=1 y 




J a=l \p=l 


y=i 


“> / "a . ,2 « 

L Il(<pJ +II 

a=l \p=l k ^ j 




ay,;' 


(5.5) 


(5.6) 


After integration over the mass we get 


oo / nl N 

a=l \p=l y=l y 


-®^(1+o2a) ‘ 
\/1 + o^A 


a=l \p=l y=l 


(5.7) 

Now to see what the total Euclidean action looks like, we move the term in square brackets into the 
exponent to give: 


(5.8) 


5total = fd^X i(3<p)'-7+ptp:p + i:;“i i(a(p+)'-7ay^ay 

+^^(l+a2A)-i + iln(l+a2A) 


with the understanding that we go back to continuum (Euclidean) space-time for the time being. 
Assuming <C nio^ and nigA < 1 (which will turn out to be equivalent to considering correlators at 
distances smaller than 1 /nio) this can be expanded in powers of 


^^^(l+o2A)-‘ + lln(l+<j2A) 


nio^A 




(5.9) 
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Note that the first term assigns a mass (equal to mo) to all the fields and fhus af fhis sfage, we have a 
theory of infinite (replicated) massive fields, where non local inferacfions (coupled fhrough powers of 
o^/ml) will be freafed perfurbafively. Then, fhe free field operafor for all the fields becomes □ + 
and all fhe inferacfion terms can be recast by functional differentiation with respect to the sources as 
usual^^: 


oo 


a=l 





+ II 

Y=1 j 



(5.10) 


Next, retaining only the leading order interaction term (corresponding to o^/ml) in ‘71/’“ gives 




1 - 


2m'i 


2 A 

mA -— 


n U^imo)[J+p] U^{mo)[J,,] 
«=i Vp=i / Vy=i 


(5.11) 


After hitting with the operators A, and then collecting all similar terms, we can finally drop the 
second index introduced for the sources i.e. set = ...=/+; = ... = 7^, to obtain 




1- 




rd 4 ) zyri 






+ 


G^mo^ 






AJ;] 


+ «-i) 


8 




z\ji 


+ L^j,v 



1 - 1 

N 

(N 

|tO 

1_1 

E,(4)'zWl 

Z[7+]Z[7+] 




+ {+ -^ —} + 2 Y^fA, 


Ed 4) Z[A1 


Ed 4 I ZIV 


'Af,v 


Z[7+]Z[7v 


n“ iZ<[7t]Z'’M7N, 


(5.12) 


where we have suppressed the mg dependence of partition functions and {+ —} represents the 
same terms as previous ones, but with time reversed (conjugated) objects Now it is apparent from 
the above equation (5.12) that corrections (of order 0{g^/ mA) to averaged free in-out correlators 
originate from and On the other hand, all of the possible partitions of the set of 

coordinates (as discussed in Sec. 4) contribute non-trivially to averaged in-in correlators. Next we 

*^It is worth mentioning that the effective action for the scalar field would have nig as mass, which is clear from the 
construction of W as in the limit a —>■ 0, the distribution in mass reduces to a delta distribution centered at nig which 
assigns mass equal to to all replicated fields with no interactions. 

'‘^Note that in the Euclidean version, there is no difference between time and time reversed objects i.e. ’+’ and 
partition functionals. But one should be careful as they are Wick (and reverse Wick) rotated differently to give different 
Lorentzian terms and hence we keep the subscripts differentiating the two types. 
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calculate the averaged in-out two point correlator and therefore only care about and 




‘W^„>lnZ[7i]-^£ ^F,kk + 




2 2 

O^nin 


+ ^I 2 


Kj 




^F,kj 


■ N 

'^^F.JlJld 


+ {^F,kj )^); 


‘W^„%^«.lnZ[7i]lnZ[72] + 


2 2 

O^ruo 


N 


lnZ[7i]£{2 

Fj [ 

\ 2 ' 


'^^F,klJld 


^F,kj 


N 




+ ^F,kj^F,kj 


-h {1 O- 2} -I- ^ ^F,kk + I ^^F,kjJl, 2 L + ^^F,kjJ2,_ 


(5.13) 


As before, the subscripts for the free propagator indicate lattice points (i.e. Ap^j^j = Apik — j)). 
Note that since A is quadratic, all odd point differentiations, i.e. differentiating an odd number of 
times with respect to sources, yield zero. Now differentiating with respect to sources, the averaged 
in-out two point correlator turns out to be 


((P(v2)(P(xi))£ ~ + KUm = ( ^2 ( ;ci X2 

+Aa^mo^ ( z z' ] (5.14) 

after taking the continuum limit of space-time. Note that only contributes to this averaged 
correlator. The subscript E again stands for Euclidean and the diagrams are to be interpreted as usual 
with the identification of solid lines with massive (nio) Euclidean propagators Ap, dots with internal 
Euclidean space-time points that are to be integrated out, and arrows specifying with-time flow i.e. 
with usual propagators. The correction diagrams are standard integrals in Euclidean QET [3, 21] and 
are given as 


i ^ 

X\ > . > X2 



1 


E (2^) 

1 


d^k 


(t2 + „2)2 = 8j2*^«K'') 


( 271)1 


ik^‘x^ 

d‘^k- 


'{B + mD 


32K^m, 


-Ki {nior) 


(5.15) 


where r is the Euclidean distance given as sj (X 2 — Ti)^ -|- {x 2 —x\Y + (y 2 —y\Y + {Z 2 —Zi)^ and 
K„{x) is the modified Bessel function of order n. Therefore in terms of Bessel functions, the av¬ 
eraged two point in-out correlator is approximately equal to 


((P(X2)(P(X1))£ 


1 

4 jj; 2^2 


niorKi {nior) 


2ml 


[mlr^Ko{mor) - mlr^Ki {nior)) 


(5.16) 


In the appendix A, we calculate the averaged in-out two-point Euclidean correlator by directly aver¬ 
aging the free massive Euclidean propagator over the distribution (5.3) and show that the two calcu¬ 
lations (leading to Eqs. (5.16) and (A.9) ) agree with each other to (at least) the first order in 


- 16 - 

































Although this closed calculation can be done for random masses since there exists a closed form ex¬ 
pression for in-out correlators for a free massive theory, it is not possible to do so in general. For in¬ 
stance if one is interested in having a distribution in the coupling constant, say P{X) Xexp (— 
for X E [0;oo), such that the standard perturbative analysis cannot be done for X > 1, the direct av¬ 
eraging of the correlators is not possible. On the other hand for g < 1 a perturbative expansion of 
our generating functionals can still be used for calculating (in general non-local) corrections to the 
averaged correlators in powers of g. However, to demonstrate that the above prescription works, 
we worked with only random masses for which a direct averaging of the correlator can be done to 
confirm its validity (See Appendix A for details). 

At last, we find that the averaged in-in one-point correlator is zero as opposed to the case in Sec. 
3.2. To see this, we only need which (upto order is equal to 


- 

nj n, 


:ln(Z[7+])ln(Z[7r])T 


2 2 


k.J 


~X^F,ki^FJcj r + {+ ^ I ^FM + ( ^.^F.kiJti 



1 

^F,kj 

F-\ 



) 

X ( —^F,kk + 




(5.17) 


and thus 


(cp(j))*(cpW)^ = 


nf nj /j /j 


= 0 


(5.18) 


Although we do not show any calculations, for the averaged in-in two point correlator (4.11) we need 
all four (2 x 2) partitions i.e. TF'V TF'V - TF'V + - and TF'V + - A functional differentiation 
of which (with respect to sources at different points as discussed) would give the desired object. 


6 Discussion 

Quantum field theories with random potentials are of importance for condensed matter systems with 
the so-called quenched disorder [1-4] and for cosmological systems in the context of string theory 
[5-7]. In this paper we made the first step towards describing such systems using the method of 
generating functionals. More precisely, we constructed two different generated functionals for cor¬ 
relators in the case of a single quantum scalar field with random potentials using the replica trick. 
The first generating functional requires a single (double) set of replicated fields only, but involves 
taking an unphysical extrapolation to negative one number of fields (See Sec. 2). The second gen¬ 
erating functional has an infinite number of fields, but involves only an interpolation between zero 
and one number of fields in each replicated set (see Sec. 4). The examples that we studied in detail 
were a random distribution of masses (see Sec. 5) and two examples of random interactions, one 
described by an (Euclidean) action of the harmonic oscillator (see Sec. 3.1), and the other described 
by a Gaussian random potential (see Sec. 3.2). We have described how the perturbation theory for 
calculating in-out and in-in correlators works for these examples; the generalization to more fields 
and to different kinds of fields should be straightforward. We have also identified a technical problem 
for the distribution of random potential described by Euclidean field theories which limits our ability 
to perform perturbative calculations analytically, but can be dealt with numerically (see Sec. 3.1). 

Before we conclude let us briefly note that the probability of a particular type of scattering event 
is the square of the corresponding scattering matrix element which is an observable quantity. Eor ex¬ 
ample, considering a one to one scattering process for a real scalar field, the corresponding scattering 
amplitude {k\p) is related to the two point correlator ((p(x 2 )tp(xi)) through the ESZ formula [20, 21], 
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but the probability is given by | {k\p) p. Then the quenched averaged amplitude of such a scattering 
process is {k\p), but the averaged probability is | {k\p) p. For systems with random potentials, one 
can argue that we must not take quenched average over the correlators/Green’s functions as they cor¬ 
respond to amplitudes for different potentials and taking averages as in (1.12) would correspond to 
letting quantum trajectories on different potentials interfere which is unphysical. Thus, one should 
rather calculate the average of probabilities of different processes as in (1.13). Although this is rea¬ 
sonable for condensed matter systems with quenched disorder, in quantum cosmology one might still 
argue that it makes sense for the wave function in different locations on the random landscape to 
interfere. To keep the discussion as general as possible we described the formalisms for calculating 
quenched averages of both amplitudes (1.12) and probabilities (1.13). Application of the developed 
methods to inflationary cosmology will be discussed in a separate publication. [12] 
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A Two-point correlator for a random mass 


Here we calculate the averaged two point in-out Euclidean correlator for random masses by directly 
averaging over the two point Euclidean correlator (Wick rotated Eeynman propagator) for a massive 
free real scalar field which is given as 
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With distribution (5.3) and working with the integral form for the propagator, we have 
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Completing squares in the gaussian integral for mass and then integrating it out gives 
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where = mj-r^ jA and 


jl. Change of variables + z = x leads to the following integral 




b^ 

dxy 1- e 


x-d 


Jb^ 

which, when expanded in powers of b^ and retaining only first order terms in square root, gives 
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This can be recast as 
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Now for b^ jcd <^\ (which is equivalent to a^/m^ <C 1) the order of the second integral is bounded 
by an exponentially small number while the first term remains finite, which leaves us 

with the first integral. Moreover we assume that b^< 1 (which is equivalent to considering distances 
r < 1 /nio) giving us 
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to first order in b^. This can be easily re-written as 
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which reduces to 
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and is the same as eq. (5.16). 
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